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I. 



INTRODUCTION 



There is an abundant literature on Coulomb systems, in particular, ions made of three 
unit charges, (mf , m^ , m^) . For the ground state with total angular momentum and parity 
L p = + , some particular configurations have been studied in great detail, and the map of 
the stability domain has been drawn, when the masses m^ are varied. For a review, and 
references, see, e.g., pp. 

The study of this + ground state indicates a number of interesting features. There is a 
variety of situations: deep binding as for H 2 + , instability as for (p,p, e~), and states at the 
edge of binding, such as Ps~(e + , e~, e~) or B.~ (p,e~ ,e~). The stability is reinforced when 
two particles become identical (m 2 = m 3 ). Schematically, the system has two components: 
a (mf, 777.3 ) cluster with a loosely attached third particle, or (mf, mj) cluster weakly linked 
to the second particle; when the thresholds become identical, the two components have 
maximal interference to build a compound system. 

For a fragile system such as H~, the Hartree-Fock picture fails. No variational wave 
function f(ri2)f(ris) gives an expectation value lower than the threshold energy. One needs 
either a explicit ^-dependence in the wave function, or a breaking of factorization, as in 
the famous Chandrasekhar wave function [2], 



in which the particle identity is first broken and then restored by a counterterm, a strategy 
sometimes named "unrestricted Hartree-Fock" j3]. In natural units, the H(ls)+e~ threshold 
for H~ is at —0.5, the above wave-function gives an energy of about —0.513, and the elaborate 
calculations about —0.528. So it is not surprising that there is only one bound state below 
this lowest threshold, as rigorously demonstrated by Hill [3]. 

However, the attention was later paid on the positive-parity state where each electron is 
a p-wave, coupled to form a total angular momentum L — 1. Such a state with L p = 1 + 
cannot decay into H(ls) + e~ as long as radiative corrections and spin effects are neglected. 
Its effective threshold consists of H(2p) + e~, at —0.125. Variational calculations, first at 



\1> = exp(-a ri2 - b r 13 ) + exp(-6 r u - a ri 3 ) , 



(1.1) 



* Electronic address: 



jean-marc.richard@lpsc.in2p3.fr 



2 




FIG. 1: Representation of a 

(m^",m^,mjj~) Coulomb system 
from its normalized inverse masses 
a, = m^ 1 /(m^ 1 + 1 + rrig x ). Some 
special cases are also shown. 



Oslo [5] and confirmed elsewhere (see, e.g., [SI [7]), give an energy of about —0.12535, just 
below the threshold. Grosse and Pittner [S] have shown that this is the only unnatural-parity 
state of H~ with this type of stability. But higher configurations of unnatural parity exist, 
for instance by attaching a positron to a two-electron system of unnatural parity [9]. 

Starting from and exchanging the masses leads to H 2 + , whose ground state is deeply 
bound and excitation spectrum very rich, and not surprisingly, includes a stable 1 + state. 
However, Mills pjj] investigated the case of equal masses, i.e., Ps~, and found the 1 + state 
to be unbound. More precisely, he studied the 1 + state for any configuration (M + , m~ ,m~) 
and got stability everywhere except for 



0.059 < m/M < 2.36 



;i.2) 



This estimate was confirmed by Bhatia and Drachman [TT]. It is the aim of the present 
note to check the domain of stability for (M + ,m~,m~) found in [TOT [TT] (without trying 
to challenge their accuracy), and to extend the study to the case of unequal masses for the 
negative charges. It is organized as follows. In Sec. [TTJ we summarize the rigorous results on 
the stability domain. The variational method and the results are presented in Sec. Ill, and 
Sec. [TVlis devoted to some conclusions. 



II. PROPERTIES OF THE STABILITY DOMAIN 

We consider the Hamiltonian 

P? Do Do 1 1 1 , 

H = -^- + -^- + -^ + — , 2.1 

2mi 2m 2 2m 3 r 12 r 23 r 23 

and focus on the lowest L p = 1 + state. By scaling, each charge can be set to unity, and one 
can impose that the inverse masses aij = 1/mj obey ol\ + ct 2 + a 3 = 1. Thus, as done for 
the + ground state [JJ H2] , each physical system with m 8 > can be represented as a point 
inside an equilateral triangle, and the inverse masses on are proportional to the distances to 
the sides. See Fig. [TJ 

In this representation, the frontier between stability and instability has the following 
properties: 

• It is symmetric with respect to the vertical axis, due to 2 <-» 3 exchange, 

• In each of the half-triangles limited by the symmetry axis, the instability domain is 
convex. 
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• The instability domain including, say, the point A 2 with inverse masses (0,1,0) is 
star-shaped with respect to A 2 . (A semi-straight line starting from A 2 enters at the 
most once a stability domain until it reaches the symmetry axis.) 

• If an energy 

E = E th (l + e), (2.2) 

is found for (M + ,m~,m~), where E t ^ = — Mm/[2(M + m)] is the (M + ,m~) thresh- 
old, and e measures the relative fraction of extra binding, and if one considers an 
asymmetric combination (m^m^m^") with mi = M and the same average inverse 
mass in the negative sector, i.e., m^ 1 + m 2 1 = 2m _1 , then stability extends at least 
up to 

l"^ 1 ~ m 2 l \ = { M ~ l +m- 1 ) . (2.3) 

The proof is the same as for the + ground state [TJ [12]. However, in this latter case, its was 
also demonstrated that any configuration with m 2 = m 3 has at most one bound state jl], so 
that the entire stability domain was connected, and clustered near the symmetry axis. This 
is not the case of L p = 1 + , as no stable state exists for equal masses [TO] [UJ. Hence, the 
stability domain includes two separate islands, one around H 2 + , and another around H~. 



III. VARIATIONAL METHOD AND RESULTS 
A. Trial wave function 



For scalar states, a generalization of (1.1) is 



^ = y^7t exp(— ax — by — cz) , (3.1) 

i 

where x = r 2 3, ...join the particles and x = |x|, ...are the relative distances. After 
integration over the trivial angular variables, one is left with integrating over xy z dxdydz, 
submitted to the triangular inequality. The matrix elements can be deduced from the generic 
function 

F(a,/?, 7 )= /// ew(-ax-0y-iz)dxdydz= {a + w A +fr)h + a) , (3.2) 

\x—y\<z<x+y 

and its derivatives. 

For a L p = 1 + state with projection L3 = j, a trial function is the superposition 

^ = (y x z)j ^ li exp(-ai x -hy - Ciz) , (3.3) 

i 

with the possibility of enforcing explicitly the x y symmetry if m 2 = m 3 . Once the integra- 
tion is carried out over the angular variables, one can express, again, all the matrix elements 
in terms of the derivatives of F . For a given set of range parameters, the weights 7^ are 
given by a generalized eigenvalue problem. It is possible to assign the cij , hi , q to have values 
in a restricted set (u, v , . . .), i.e., (a^, bi, q) = (u, u, u), (u, u,v), . . . (u,v,w), . . . (w, w, w), . . ., 
and even assume (u, v, . . .) in a kind of progression, in order to simplify the minimization 
over the non-linear parameters, without a significant loss of accuracy. This is similar to the 
strategy used, e.g., by Kamimura et al. when handling the expansions over Gaussians [IB"] . 
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B. H -like states 



For H~, the energy is found to be E ~ —0.1253, as in the literature. For other symmetric 
configurations (M + ,m~ ,m~), we obtain stability for M/m > 2.39, in good agreement with 
Mills's estimate (1.2). The width of the stability domain is very narrow. For an infinitely 
massive proton, we find stability only for m^/m2 ^5 1.006 (if m% > rri2). This is very close to 
the estimate based on (2.3), as expected for such a small width. For instance, the 1 + state 
does not exist for the very exotic (p, tt~ , fi~) system. 



C. H9 + -like states 



Stability is found for (M + ,m ,m~ ") with M/m < 0.055, again in good agreement with 
Eq. ( 1.2 ). On the side of the triangles, we obtain stability for the inverse masses (1 — a, a, 0) 
and the symmetric points (1 — a, 0, a) for < a ^ 0.094. For the point of the symmetry 
axis with the inverse masses (1 — a, a/2, a/3), the relative excess of energy is found to be 



0.034, and Eq. (2.3) is underestimates the width of the domain by about 15%. 




FIG. 2: Domain of stability for (m^m^mj) in the triangle of normalized inverse masses: the 
ground state + (dashed lines) and the unnatural-parity state 1 + (solid). On the right, magnifica- 
tion of the bottom part of the plot. 



D. Summary 

The domain of stability contains at least the areas displayed in Fig. [2] The frontier is 
almost straight, so the convexity of the instability domain is weakly pronounced. The spike 
around H~ is very small and extremely narrow, and shows up only after magnification. For 
comparison, the domain of stability of the ground state + is also shown. 
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IV. OUTLOOK 

The case of unnatural-parity states of three unit-charge particles displays an interesting - 
and rare - example of discontinuous stability domain, where stability disappears and shows 
up again, when some constituent masses are varied continuously and monotonously. This 
means that the transition from a molecular type of binding (H2 + ) to a halo-type of binding 
(H~) involves some more fragile intermediate dynamics of stability. 
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